ABSTRACT 


A consistent  environmental-acoustic  model  for  a deep  moving 
ocean  is  formulated.  The  acoustic  model  for  regularly-perturbed 
SOFAR  rays  is  approximately  solved  using  a type  of  WKB(J)  expansion. 
Interfacing  conditions  between  the  hydrodynamics  and  acoustics 
are  developed  which  lead  to  constraints  on  acoustic  frequency  and 
transmission  range.  As  an  application,  transmissions  are  con- 
sidered through  stochastic  internal-wave  fields,  which  have  been 
modeled  in  a previously-published  paper  by  the  authors.  Formulas 
for  ray  phase  variances  are  derived.  These  formulas  are  asymptot- 
ically evaluated  for  rays  with  relatively  significant  depth 
variation,  using  the  stationary  phase  method.  New  results  are 
obtained  for  the  dependence  of  the  variances  on  internal-wave 
primitives,  such  as  energy  spectra.  Expected  multipath  intensity 
is  calculated  for  transmission  through  an  ocean  with  static  state 
modeled  by  a bilinear  sound-speed  profile.  The  effects  of  the 
internal-wave  field  and  of  varying  internal-wave  parameters  on 
the  expected  intensity  are  shown  to  be  significant. 


INTRODUCTION 


In  recent  years,  oceanic  environmental-acoustic  models  have 
been  developed  to  consider  a variety  of  phenomena,  such  as  currents 
and  internal  waves.  The  objective  of  most  of  these  studies  is  to 
predict  variations  in  acoustical  quantities,  as  for  example  intensity, 
when  sound  transmissions  traverse  an  ocean  region  containing  one  or 
more  environmental  effects.  The  first  principal  purpose  of  this 
paper  is  to  formulate  a consistent  environmental-acoustic  model 
incorporating  both  sound-speed  and  velocity  variations.  Primary 
features  of  the  model  include  its  orientation  toward  deep-ocean 
transmissions,  its  foundation  in  ray  acoustics,  and  its  suitability 
for  both  deterministic  and  stochastic  sound-speed  and  velocity 
fluctuations.  A central  and  unique  characteristic  of  the  model  is 
its  emphasis  on  consistency.  By  a consistent  analytic  model,  we 
mean  one  for  which  the  model  equations  consistently  include  all 
terms  and  effects  of  a specified  asymptotic  order,  and  exclude  all 
others  of  smaller  order.  Consistency  must  be  observed  in  the 
environmental  and  acoustical  portions  of  the  model  separately  and 
in  their  relationship,  or  interfacing.  As  one  example  of  the 
consequences  of  consistent  interfacing,  estimates  can  be  found 
for  the  applicability  of  acoustic  predictions,  based  on  approximations 
made  in  determining  expressions  for  sound-speed  and  velocity. 

Our  second  principal  objective  is  to  apply  our  model  to 

investigate  features  of  sound  transmission  through  a stochastic 

field  of  internal  waves.  Although  much  recent  work  has  focused 

on  this  problem,  many  aspects  have  been  incompletely  resolved  or 

overlooked  entirely.  For  example,  we  know  of  only  two  acoustical 
1 2 

treatments  ' which  include  velocity  as  well  as  sound-speed 
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fluctuations;  the  effects  of  the  former,  and  the  circumstances  in 

which  they  are  important,  have  not  been  sufficiently  explored. 

A second  example  is  our  investigation  of  internal-wave  effects 

on  expected  multipath  intensity,  which  is  evidently  unique  in  the 
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literature.  As  a third  example,  most  analytic  investigations 
have  used  an  internal-wave  model  first  introduced  in  Ref.  10  and 
updated  subsequently.  The  use  of  inconsistent  approximations  in 
that  model,  and  acoustic  work  based  on  it,  can  be  shown  to  signifi- 
cantly alter  predictions  for  the  acoustic  quantities  we  investigate 
here.  We  note  that  our  consistent  approach  differs  from  acoustic- 
fluctuation  models  that  neglect  explicit  coupling  to  environmental 
effects  (for  example,  as  in  portions  of  Refs.  11  and  12).  Also, 
other  acoustic  internal-wave  studies13'14  have  been  directed  toward 
parameter  regions  where  ray  theory  may  not  be  appropriate.  At 
the  beginning  of  Sec.  Ill,  we  have  listed  six  points,  on  at  least 
one  of  which  our  development  differs  from  each  previous  investigation. 

In  Sec.  I,  a ray  acoustic  model  for  moving  media  is  formulated. 
Equations  of  the  model  are  solved  using  a WKB(J)-type  approximation 

under  a number  of  explicitly  stated  assumptions.  The  application 
of  the  model  to  both  deterministic  and  stochastic  current  and 
sound-speed  fluctuations  is  discussed.  Appropriate  variable  scales 
and  model  parameters  are  selected  in  Sec.  II.  Consistency  and  inter- 
facing conditions  on  the  acoustic  and  hydrodynamic  parameters  are 
discussed,  and  their  constraints  on  CW  source  frequency  and  trans- 
mission range  are  illustrated.  Useful  simplifications  based  on 
consistency  are  presented.  In  Sec.  Ill,  we  derive  formulas  for 
ray  phase  variances  in  the  presence  of  a stochastic  internal-wave 
field.  Results  are  heavily  based 
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on  the  consistent  internal-wave  model  presented  by  the  authors  in 

15 


a previous  paper.  We  define  a class  of  rays  with  relatively 
significant  depth  variation,  and  asymptotically  evaluate  the 
phase  variance  for  such  rays  in  Sec.  IV  using  the  stationary 
phase  method.  Prom  these  results,  we  extract  the  specific  internal- 
wave  parameters  which  influence  phase  variance,  discuss  their 
effects,  and  compare  our  predictions  to  previous  calculations. 

In  Sec.  V,  we  apply  our  phase  variance  results  to  the  calculation 
of  multipath  expected  intensity,  assuming  a bilinear  static 
sound-speed  profile.  The  effects  of  the  internal-wave  field  and 
of  varying  internal-wave  parameters  on  the  range  dependence  of 
intensity  are  illustrated.  Finally,  our  results  are  summarized 
in  Sec.  VI. 
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I.  ACOUSTIC  FORMULATION  AND  OCEAN  APPLICATION 

We  consider  the  acoustic  pressure  p(x,t),  as  a function 
of  nondimensional  position  x and  time  t,  in  a medium  with 
nondimensional  sound  speed  c(x)  and  velocity, u (x)  but  with  no 
attenuation  or  dispersion  effects.  Assuming  the  medium  is  unaffected 


by  acoustic  signal  propagation,  the  nondimensional  equation  for  p is 

I 2 
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8/at  + M(u*7)j’  p = c2  V2p  , 


(1) 


where  7 (V  ) is  the  gradient  operator 
(Laplacian) . The  dimensional  position, 

sound  speed,  velocity,  and  time  are  nondimensionalized  by  the 
characteristic  values  XQ,  cQ,  UQ,  and  XQ/c0*  respectively,  all 
to  be  specified  later.  The  Mach  number  M = Uq/°q  is  the  ratio 
of  the  characteristic  velocity  to  the  characteristic  sound 
speed  and  is  assumed  small.  It  is  appropriate  to  expand  the 
total  sound  speed  c as  a sum  of  a static  term  c and  a velocity- 
induced  term  Me,  which  we  further  expand  in  M: 


c + M[c*°*  + Me*1*  + ...] 


Similarly,  the  velocity  u has  the  expansion 


-*■  -WO)  -*•  ( 1 ) 

u=uv  ’ + Mu ' ' + ... 


(2a) 


(2b) 


The  time  dependence  is  separated  from  the  spatial  dependence 
in  Eq.  (1)  by  assuming  p to  be  of  the  form 


p(x,t)  * Re  (x)  exp(-iu)t)]  , 


(3a) 


where  Re  denotes  real  part,  i is  the  imaginary  unit,  and  o>  is 


U 
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dimensionless  angular  frequency* 

Substituting  Eq.  (3a)  into  Eq.  (1)  , we  obtain  the  reduced  wave 
equation  for 

-2  2 -*■  2 

c V \|»  + [o)+iM(u*7)]  = 0 . (3b) 

We  restrict  our  attention  to  large  values  of  u>,  the 
geometric  acoustics  limit.  We  assume  an  asymptotic  relation 
between  M and  w \ M = u>  ^G(u  ^)  , such  that  M is  asymptotically 
order  one  with  respect  to  u>  1;  i.e.  as  u"1  + 0,  g + F where  F 
is  a bounded  constant.  This  assumption  has  been  used  for  other 
acoustics  problems,  as  in  Ref.  17  for  example.  Substituting 


into  Eqs.  (2)  and  (3b) , the  Mach  number  expansions  of  sound 
speed  and  velocity  reduce  to  expansions  in  oj.  To  obtain  an 
approximation  for  ip,  we  use  a type  of  WKB(J)  expansion18  in  the 
parameter  u 1.  We  assume  that  the  complex  function  ip  has  the 
polar  form 


i|<(x)  = a(x)  exp  [i4>  (x)  ] , 


(4a) 


and  expand  the  real  functions  amplitude  a(x)  and  phase  4>(x)  as 
follows : 


(0) 


-1  (1) 
o)  a 


(4b) 


and 

* = U)T  + + U)_1*(1)+  ...  . (4c) 

The  function  ux  (x)  is  the  principal  contributor  to  the  phase  of  ij>. 
The  terms  u ^a^and  u>  ^4^  are  order  oj  ^ corrections,  and  these 
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provide  approximations  to  ^ and  p correct  to  order  w”3 . Equations 
for  x , a^,  and  4>  ^ are  obtained  by  substituting  Eqs.  (2)  and 
(4)  into  Eq.  (3b)  and  equating  both  real  and  imaginary  parts  of 
like  powers  of  to  to  zero. 

We  now  focus  on  the  order  one  solution  for  p,  with  brief 
consideration  of  higher  order  terms  deferred  to  Sec.  II.  For 
convenience,  variables  c,  u,  Q,  and  a without  superscripts 
hereafter  represent  those  variables  previously  written  with 
superscript  zero.  The  equations  obtained  for  t,  a,  and  <J>  are 

|Vt|  = c”1  , (5a) 

s 

V»(a2Vx)  = 0 , (5b) 

and 

O "I 

Vt*V4  + F c I c c + u*Vt]  = 0 . (5c) 

s s 

The  surfaces  of  constant  phase  are  t(x)  = constant,  and  curves 
normal  to  these  surfaces  are  the  rays.  Each  ray  describes  a 
curve  in  space  as  a parameter,  which  we  choose  as  arclength  s, 
varies.  Moreover,  each  ray  is  a member  of  a two-parameter  family 
of  curves,  which  we  emphasize  by  denoting  a ray  as  x = x(s;£,n)« 

The  variables  (;  and  n parameterize  the  initial  constant-phase 
surface,  i.e.  s = 0,  from  which  all  rays  emanate.  In  this  paper, 
we  consider  only  rays  that  are  continuously  refracted,  and  mention 
later  the  modifications  that  arise  if  boundary  reflections  occur. 

In  addition,  we  note  the  implicit  assumption  that  the  presence 
of  the  order  M sound  speed  and  currents  produces  a perturbation 
in  a ray  of  regular,  as  opposed  to  singular,  type.  For  example, 
we  forbid  the  small  fluctuations  to  change  a ray  from  continuously 


: 


MMMMi 


I 

refracted  to  refracted/boundary-reflected  type,  or  to  cause  a 

it 

ray  to  disappear  altogether. 

..  Once  the  equation  of  a ray  is  known,  the  traditional  form  of 

the  ray  travel  time  T (s; e , n) =t [x (s; £ , n) ] is  obtained: 


In  eq.  (6a),  the  arguments  of  the  integrand  are  l,  £,  and  n. 
Similarly,  we  obtain  the  order  one  phase  along  a ray, 
♦(s;E,n)  = *Ix(s;C,n)]: 


i: 

i: 

[ 

o 

n 


- F 


r 


(c  2 (c  + e*u) }dt  , 

S 


(6b) 


where  <^(s;£,n)  = 3x(s;£,n)/3s  is  the  unit  tangent  vector  to  a ray. 

In  order  to  find  an  expression  for  amplitude  along  a known 
ray,  A(s;£,n)  * a [x(s; £ , n) ] , we  consider  the  important  example  of 
a point  source.  One  model  for  a point  source  is  a spherical 
initial  phase  surface  of  nondimensional  radius  r = 1.0  m/XQ. 

The  pressure  on  this  sphere  is  assumed  a constant  value  which  we 
choose  as  the  scale  for  the  nondimensional  pressure,  so  that 
A(0;C,n)  = 1 • The  parameters  5 and  n are  respectively  taken  as 
latitude  and  longitude  on  the  sphere,  so  that  -ir/2  <_  C <_  n/2 
and  0 <_  n <.  2ir.  For  continuously-refracted  rays  emanating  from 
a point  source  centered  at  x = x^  , the  amplitude  can  be  shown 
to  be 


A ■ r[cs(cos  S)/c^  J] 


1/2 


(7a) 


X. 


-8- 


In  Eq.  (7a)  c ^ = cg(x^  )»  and  geometrical  spreading  loss  is 

incorporated  through  J,  the  magnitude  of  the  Jacobian  of  the 
ray  parameterization, 

. -►  . "f 

, 9x  9x  -*■  . 

J ■ 35  x 57;  • e • (7b) 

The  order  one  total  acoustic  pressure  at  a receiver 
located  at  x = x is  constructed  next,  assuming  a fixed  CW  point 

<R. 

source  at  frequency  u>.  If  P(s,t;^,n)  = p [x(s;  £,n)  ,t]  is  the 
order  one  acoustic  pressure  along  a ray,  then  from  Eqs.  (3a) 
and  (4) 

p = A cos  {w(T-t)  + (j)}  . (8a) 

The  order  one  total  pressure  at  the  receiver,  p (x^  ,t)  is  the  sum 
of  the  contributions  from  individual  rays.  Let  the  set  A be 
defined  as 

A = { (s,£»n)  |x(s;£,n)  = xc  } , (8b) 


and  assume  members  of  A can  be  labeled  by  integer  values  of  an 
index  n.  If  P(nj  is  pressure  corresponding  to  the  n th 
member  of  A,  then 


(Xp/t)  ]CP(n)  * 


Other  acoustic  quantities  of  interest  follow  directly  from 
Eqs.  (8).  For  example,  the  total  amplitude  a(Xp)  and  total 
phase  5(3  ) at  the  receiver  are  related  to  p by  Eqs.  (3a)  and 
(4a)  : 


p(x^,t)  = Re{a(x^)  exp[i(5(x^)  — a>t)  ] > 


* 


Thus  a , which  is  proportional  to  the  acoustic  intensity,  is 


-2 

a 


£A(n>  <*Pli<»T(n)  + *(n))I 


(9b) 


To  specialize  our  moving-media  acoustic  formulation  to 
deep-water  oceanic  transmission,  we  choose  a Cartesian  coordinate 
system  with  origin  on  the  ocean  surface.  The  positive  x,y,  and  z 
axes  are  in  the  directions  of  east,  north,  and  increasing  depth, 
respectively,  with  unit  vectors  t,  j,  and  £ (Fig.  1) . Also 
shown  in  Fig.  1 are  the  location  of  the  source  ^ at  x^  = (x ^ ,y^  ,z^ ) 
and  the  parameters  £ and  n*  The  latter  are  spherical 
coordinates  with  zero  values  in  the  plane  z = z^  and  y = y^  , 
respectively,  and  with  increasing  £ values  in  the  direction  of 
increasing  z. 

The  static  sound  speed  in  the  ocean  is  well  approximated  by 


a function  of  depth  only,  c = c (z).  This  simplification  leads 

s s 


to  the  first  integral  of  the  ray  equations, known  as  Snell's  Law: 


cos  6 / cos  i = cg(z)  / Cj  , (10a) 

where  9 is  the  inclination  angle  of  the  ray,  measured  in  the 

same  sense  as  £ (see  Fig.  1).  Since  we  consider  only  continuously 

refracted  rays  between  fixed  source  and  receiver,  it  can  be  shown 

that  a unique  value  n = n exists;  i.e.  rays  are  confined  to  the 

s 

vertical  plane  containing  source  and  receiver.  We  accentuate 

this  symmetry  by  eliminating  arclength  as  the  ray  parameter  in  favor 

2 2 2 

of  range  R from  the  source,  R = (x  - ) + (y  - ) . If  we 

regard  z = z(R ;£)  and  0 = 0(r;£)  along  a ray,  it  follows  that 
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(dz/dR)2  = tan2e  = [c£(sec£  )/Cg  (z)  ] 2 - 1, 

z(0;U  = z^|  . (10b) 

We  note  that  the  set  A from  Eq.  (8b)  may  be  expressed  as 


A = { (C»ns)  |z(x»C)  “ z(l  • tan  hs  = 

<y«.'  yi  )/<x«-"  XJ,)  ’ 


(10c) 


where  the  dimensionless  parameter  x equals  the  dimensional 
horizontal  range  RQ  from  source  to  receiver  divided  by  XQ. 

We  obtain  the  leading  approximation  to  the  travel  time  of 
a ray  to  the  receiver  from  Eq.  (6a) , after  changing  integration 
variable: 

X 


T(X,U  = c.  (sec  5)1  c 2 dR 


1 


(Ha) 


From  Eq.  (6b) , we  find  similarly  that  the  order  one  phase  at 
the  receiver  is 

X 

4>(X?£rn)  = - F (sec 5)  J c~3  (c  + e*u)dR,  (lib) 

0 

where  the  unit  tangent  vector  to  a ray  is 

e = (cos9)  [ (cosn  )i  + (sinnc)D]  + (sine)lt  . (He) 


Note  that  <J>  is  independent  of  n_if  the  functions  c and  u are 

symmetric  about  the  z-axis.  Both  integrals  in  Eqs.  (11a)  and 

(lib)  are  evaluated  using  the  ray  paths  z = z(R;£)  as  solutions 

of  Eqs.  (10).  The  function  J in  Eq.  (7b)  is  invariant  under 

19 

the  change  of  parameter  from  s to  R.  It  follows  from  Eqs.  (7), 
(10)  , and  (11c)  that  the  amplitude  at  the  receiver  is 

-1/2 

A(x;5)  = r | x ( 3 z/3 5 ) | 


(Hd) 
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In  this  section,  we  have  assumed  that  rays  are  continuously 
refracted  and  never  reflected.  When  reflections  occur,  the  ray 
can  be  considered  a succession  of  segments  between  the  reflections. 
Within  each  segment,  our  approach  may  be  extended  readily. 

Analogous  formulas  are  found, but  with  appropriate  modifications, 
one  example  being  the  limits  of  integrals  in  Eqs.  (11) . 

However,  boundary  models  must  be  specified,  along  with  their 
effects  on  acoustic  amplitude  and  phase. 

The  foregoing  acoustic  formulation  is  useful  for  either 
deterministic  or  stochastic  sound  speeds  and  velocities.  In 
applications  involving  internal  waves,  for  example,  these  quantities 
are  often  regarded  as  members  of  statistical  ensembles.  We  have 
assumed  here  only  that  c and  u depend  on  x.  Precise  specification 
of  the  ensemble  depends  on  the  particular  application. 

From  Eq.  (lib) , these  stochastic 

processes  induce  a stochastic  process  for  the  phase  4,  which  in 
turn  makes  the  total  received  pressure  random.-  One  example  of  a 
statistic  of  the  received  pressure  is  the  expected  intensity,  with 
expectation  defined  here  as  an  ensemble  average  and  denoted  by  < > . 

A formula  for  this  quantity  for  general  multipath  propagation 
is  obtained  from  Eq.  (9b)  : 

<«2>  = EA<„>  + E EA(„)A(m)cos  I“(T(n)-T(m))1  ' 

n n m 

n/m 

“ exp  5 • <12a> 


The  expectation  on  the 
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right  side  of  Eq.  (12a)  is  equal  to  the  characteristic  function 

of  the  random  variable  6,  . - $ . . evaluated  at  one. 

(nj  (m) 

Simplification  of  Eq.  (12a)  is  possible  if  specific 

assumptions  are  made  about  statistical  properties  of  the  $'s. 

Examples  include  assuming  that  $ ^ and  <p  ^ are  independent  or 

that  d> . . and  . have  bivariate  normal  distributions.  With 
(n)  (m) 

both  of  these  assumptions,  and  using  in  addition  < 4>  > = 0, 
it  follows  that 

< 0,2  > = E *2n)  + I E A(n)A(m)COsI“(T(n)"T(m))] 
n n m 

n/m 

exp  t_(a2m)  + °2n))/21  ' (12b) 

where  we  have  expressed  the  characteristic  function  of  a 

normally-distributed  random  variable  in  terms  of  its  variance, 

2 2 

a = < <f>  >.  The  first  term  in  Eq.  (12b) 

is  proportional  to 

intensity  either  if  multipath  interactions  are  ignored,  i.e.  if 

20 

incoherent  multipath  summation  is  employed,  or  if  the  phase 
variances  are  assumed  large,  in  the  absence  of  phase  variances, 
the  second  term  reduces  to  the  contribution  of  interacting 
multipaths  to  the  intensity.  When  phase  variations  are 
present,  each  term  in  the  second  sum  is  reduced  by  a damping  factor 
which  varies  in  magnitude  depending  on  the  ray  pair  involved  in 
the  term. 


II.  CONSISTENCY  AND  SCALING 

In  this  section,  we  investigate  the  application  of  our 
acoustic  model  to  deep-water  ocean  transmission  problems.  In 
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particular,  we  are  concerned  with  interfacing  the  acoustic 
model  and  an  appropriate  hydrodynamic  model  for  ocean  current 
and  sound-speed  distributions.  This  interfacing  involves  two 
steps.  First,  appropriate  scales  XQ,  UQ,  and  cQ  are  selected, 
and  the  magnitudes  of  the  dimensionless  parameters  u and  F 
estimated.  We  then  verify  that  these  parameters  satisfy  the 
asymptotic  restrictions  imposed  by  the  acoustic  model.  A second 
interfacing  step  is  necessary  in  many  applications  (see,  for 
example.  Refs. 15  and  21)*  This  is  because  the  sound  speed  and 
current  are  obtained  from  asymptotic  approximations  of  the  hydro- 
dynamic  equations. 

We  first  briefly  consider  ocean  thermodynamics,  in  order 
to  introduce  notation  and  to  obtain  a consistent  and  general 
formulation  for  the  static  state.  Sound  speed,  density# and  pressure 
are  thermodynamic  variables  related  by  an  equation  of  state,  which 
follows  from  relations  in  Ref.  22,  for  example.  Such  empirical, 
equations  are  typically  analytically  cumbersome,  but  consistency 

permits  some  simplification.  Let  the  dimensional  sound  speed  c*, 

* * 

density  p , and  ocean  pressure  !r  be  related  by  the  generic 
* * * * 

state  equation  c = c (p  , p ) , assuming  constant  salinity. 

Further,  let  pQ,  IP0»  and  c0  = c*(pQ,  f*o^  be  characteristlc 

values  of  these  variables  at  the  ocean  surface,  with 

approximate  values  pQ  = 1.02  x 103  kg  m~3,  (p  » 1 atm,  and 

cQ  * 1540  m sec”1.  Analogous  to  the  definition  of  c in  Sec.  I, 

we  define  nondimensional  density  p and  pressure  (p  by  p*  = pQ  p 

and  P - PQ  = (pQgD)  <p  , where  g is  the  gravitational 

_2 

constant,  9.81  m sec  . Expanding  the  equation  of  state  in  a 


-14- 


I 


Taylor's  series  about  p*  = pQ  and  P*  = (?Q  and  nondimensionalizing 


yields 


c = 1 + v(l-p)  + Ayp  + ...  , 


(13a) 


where 


v = ■ Poco1(ac*/»P*)o 


y = P0c0  Oc*/3P  ) 


(13b) 


(13c) 


The  zeros  in  Eqs.  (13b)  and  (13c)  denote  evaluation  at  (Pq»Pq)» 

For  example,  the  values  of  v and  y obtained  from  the  state  equation 

of  Ref.  22  fora  35  o/oo  isohaline  ocean  are  5.3  and  7.9, 

respectively.  We  remark  that  although  these  numbers  may  change 

with  different  oceanic  conditions,  i.e.  other  values  of  p^,  Cq, 

and  salinity,  the  value  of  y - v remains  approximately  constant. 

2 

The  parameter  A = gD/cQ  appearing  in  Eq.  (13a)  is  the  ratio  of 

the  hydrostatic  pressure  scale  pQgD  to  the  thermodynamic  scale 
2 

pQCg  . For  a deep  ocean  (D  = 5 km) , A is  a small  parameter, 
with  approximate  value  2.1  x 10  . 

Specification  of  exactly  one  of  the  three  thermodynamic 
variables,  as  a function  of  dimensional  depth  z*,  determines 
the  others  via  the  state  equation  and  hydrostatic  balance.  One 
way  is  to  prescribe  the  static  density,  which  is  conveniently 

* It 

written  as  follows.  If  we  let  z = z /D  be  depth  nondimension- 
alized  by  total  depth  0,  then  the  nondimensional  static  density  is 


p (z)  = 1 + A(z  + ir(z) ) , 


(14a) 
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where  n(z)  is  the  potential  density.  Substituting  Eq.  (14a) 

p.  A 

and  the  hydrostatic  result  P = z + 0(A)  into  Eq.  (13a),  we 

s 

obtain  the  static  sound  speed  as 

c ( z ) = 1 + As  (z)  , (14b) 

s 

where  s to  order  one  in  A is 

A A 

s(z)  = (y-v)z  - vn(z)  . (14c) 

The  small  parameter  A in  Eqs.  (14a)  and  (14b)  accentuates  the 

weak  dependence  of  p and  c on  depth.  He  note  that  for  some 

s s 

acoustic  applications,  it  may  be  more  convenient  to  specify  the 

A 

sound-speed  variation  s(z),  from  which  Eqs.  (14a)  and  (14c) 

may  be  used  to  find  p . Equations  (14)  are  written  in  terms  of 

8 

A 

the  nondimensional  Brunt-Vaisala  frequency  N(z)  (scaled  by 
g/c_),  using  the  identity 

0 A 

z 

Tf(z)  = J N2(s)ds  . (14d) 

0 

In  fact,  for  an  exponential  N(z),  Eqs.  (14b)  - (14d)  reduce  to 

23 

a special  case  studied  earlier.  Our  general  formulation  is  appro- 
priate for  arbitrary  specification  of  the  ocean  static  state.  Also, 
it  is  readily  generalized  to  include  salinity  variations,  which 
are  not  considered  in  this  paper. 

Turning  to  specification  of  the  scales,  we  consider  the 
length  scale  XQ  first.  The  ocean  depth  D was  used  for 

A 

convenience  in  scaling  the  variable  z in  the  thermodynamic 
formulation.  However,  the  appropriate  XQ  is  actually  the  minimum 
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of  the  plausible  choices  in  the  environmental-acoustic  problem. 

The  most  obvious  of  these  possibilities  are,  besides  D,  the 

range  of  transmission  and  the  scale  height  of  the  dimensional 

* 

static  sound  speed  c . if  we  assume  hereafter  that  Rft  is  at 

s 0 

least  as  large  as  D,  then 


X0  = Min  {D,  max  |c*|/max  |dc*/dz*|} 


(15a) 


We  may  obtain  an  estimate  for  the  scale  height  which  is 
independent  of  the  particular  functional  form  of  c*.  Using 
Eqs.  (14b)  - (14d) , Eq.  (15a)  becomes  approximately 


XQ  * Min 


in  |d,  cjjg-1  ] 


Max  | (y-v) -vn  (z)  | 


n • 


(15b) 


Estimating  N from, for  example.  Ref.  24,  the  quantity  in  braces 
is  invariably  less  than  ten,  so  the  minimum  scale  height  is 
approximately  25  km.  For  the  deep-ocean  value  D * 5 km,  it 


follows  that  XQ  * D,  z = z , 


and  X * RQ/D.  s*nce 


is  prescribed  by  the  appropriate  hydrodynamic  model,  typically 
in  value  between  1.0cm  sec  *and  ]j0  m sec  all  scales  are  now 
specified.  We  remark  that  if  other  physical  processes  are 
considered,  additional  possibilities  may  arise  for  the  choice  of 
XQ.  For  example,  models  for  scattering  effects  of  random  sound- 
speed  and  current  fluctuations  lead  to  correlation  lengths  which 
may  influence  the  choice  of  XQ,  depending  on  their  magnitudes. 

We  next  consider  restrictions  on  the  dimensionless  parameters 
u and  F.  Our  asymptotic  acoustic  model  is  valid  for  large 
values  of  w.  in  terms  of  dimensional  frequency  and  wavelength 
of  the  acoustic  signal,  f ■ c u/2«d  and  X * cft/2nf  , this 
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requirement  is  equivalent  to  D >>  A^.  This  constraint  forbids 

diffraction  phenomena,  one  of  the  traditional  limitations  of 

geometrical  acoustics.  We  note  that  the  conditions  for 

geometrical  acoustics  imposed  in  Ref.  25  are  obtained  by  choosing 

Xq  equal  to  the  scale  height  of  the  static  sound  speed  and 

requiring  u>  >>  1.  Thus,  our  choice  of  smaller  XQ  insures  that 

conditions  in  Ref.  25  are  satisfied.  Having  specified  u> 

and  M = Uq/Cq,  we  choose  for  simplicity  u>M  = G(0)  = F.  With 

2 

XQ  = D,  it  follows  that  F = 2TrfgDU0/c0.  Further,  we  require 
that  the  magnitude  of  F is  approximately  unity  or  less,  since 
otherwise  M could  be  assumed  other  than  order 

We  now  define  e as  the  order  of  the  maximum  tolerable  error 
we  accept  using  our  order  one  asymptotic  approximation  for 
acoustic  pressure,  Eqs.  (8).  An  estimate  of  the  error  incurred 
in  using  Eqs.  (8)  is  the  order  of  the  first  term  neglected  by 
the  approximation;  i.e.  w-1p^  (x,t)  by  Eqs.  (4).  in  general, 
it  can  be  shown  that  p^  has  secular  behavior  with  respect  to 
the  variable  x.  Moreover,  these  secular  errors  are  of  the  order 
of  the  greatest  nondimensional  ray  arclength  between  source  and 
receiver.  Since  we  consider  only  SOFAR  rays  here,  it  follows  from 
calculations  leading  to  Eq.  (10b)  that  the  order  of  the  secular 
terms  in  p*1*  is  x*  Therefore,  the  order  of  terms  neglected  in 
our  approximation  reduces  to  ^ • Using  the  definition  of  e, 
we  have 

0(x«_1)  < 0(e)  , (16) 

and  we  refer  to  Eq.  (16)  as  the  secularity  condition,  we  note 
that  since  x > u » 1 as  previously  required. 
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Before  illustration  of  Eq.  (16),  we  turn  to  the  consistency 
between  the  acoustic  and  hydrodynamic  models.  There  are  two 
aspects  of  this  interfacing  step.  First,  we  must  ensure  that  the 
hydrodynamic  effects  contribute  to  the  order  one  acoustic 
pressure  at  an  order  greater  than  the  acoustic  error  e.  Other- 
wise, it  is  inconsistent  to  include  them  in  the  acoustic  analysis. 
Second,  the  derivation  of  the  acoustic  model  implicitly  assumes 
that  the  current  u and  motion- induced  sound  speed  c are  known  as 
order  Mach  number  quantities.  However,  the  hydrodynamic 
derivations  of  models  for  these  variables  neglect  small  terms. 

We  denote  the  order  of  the  largest  such  neglected  terms  by  y. 
Therefore,  we  must  insure  that  neglect  of  order  y terms  in  the 
hydrodynamic  model  is  consistent  with  a tolerable  error  of 
order  e in  the  order  one  acoustic  pressure. 

To  quantify  the  above  requirements,  we  note  that  the 

hydrodynamic  variables  c and  u occur  in  the  acoustic  model  only 
in  the  integrand  of  the  order  one  phase  $,  Eq.  (lib).  Estimating 
the  integral  of  Eq.  (lib)  shows  that  the  order  one  and  order  y 
hydrodynamic  effects  produce  contributions  to  $ of  order  xF 
and  order  y\F,  respectively.  Using  this  result  in  Eqs.(8) 
and  expanding  about  xF  = 0»  we  see  that  order  one  hydrodynamic 
effects  are  consistently  included  provided 


0(xF)  > 0(e) . 


(17a) 


Similarly,  an  expansion  about  yxf  * 0 yields  the  following 
inequality  for  the  consistent  neglect  of  order  y hydrodynamic 
effects: 


■amm 


— JU~.  — 
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O(YXF)  < 0(e)  . (17b) 

We  refer  to  Eqs.  (17)  as  the  consistency  conditions,  noting 
that  an  expression  for  F has  been  determined  previously. 

For  any  given  model  application,  g,  D,  cQ,  and  UQ  are 

specified.  Thus,  the  secularity  and  consistency  conditions 

place  constraints  on  acoustic  frequency  and  transmission  range. 

We  illustrate  these  first  for  the  choices  y = e = A (roughly 

one  per  cent) . As  motivation,  we  note  that  y = A is  often  a 

useful  estimate  for  hydrodynamic  derivations  (see,  for  example. 

Ref.  21  ) . indeed,  the  importance  of  the  parameter  A is 

suggested  by  Eqs.  (14).  Since  our  parameters  are  referred  to  a 

power  of  ten,  we  replace  the  order  inequalities  Eqs.  (16)  and 

26  —1 

(17)  by  the  order  of  magnitude  comparisons  < 3A, 

XF  >_  A/3,  and  xF  ± 3.  In  Fig.  2,  we  graph  these  inequalities, 
along  with  the  earlier  assumption  x > 1,  in  dimensional 
frequency-range  space.  The  previous  estimates  of  g,  d,  and  cQ 
are  used,  as  well  as  UQ  values  of  1,  3,  and  10  cm  sec-1.  For 
each  Uq,  the  simultaneous  solution  of  the  four  inequalities  is 
the  interior  of  a trapezoid,  which  gives  approximate  ranges 
of  experimental  operating  conditions.  The  UQ  - dependent 
parallel  sides  are  determined  from  the  consistency  conditions, 
and  the  non-parallel  sides  from  x * 1 and  the  secularity 
condition.  A wide  range  of  experimental  operating  values  of 
Rq  and  fp  fall  within  each  "interface  trapezoid".  For 
example,  for  * 1 cm  sec  1 and  fQ  = 100  Hz,  the  permissible 
ranges  are  approximately  5 to  650  km.  Similarly,  the  range  of 
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frequencies  for  Uq  = 3 cm  sec  ^ and  R^  = 200  km  is  about 
30  to  560  Hz.  We  note  that  several  values  of  Uq  are  shown, 
because  among  the  environmental  parameters,  this  one  has  the 
largest  possible  relative  variation.  As  shown  in  Fig.  2, 
decreasing  values  of  UQ  typically  increase  the  allowable 
frequencies,  for  fixed  range.  This  results  because  smaller 
velocity  scales,  i.e.  smaller  Mach  numbers,  are  associated  with 
smaller  fluctuations  that  can  be  resolved  by  higher  acoustic 
frequencies. 

To  illustrate  effects  of  changing  the  tolerable  acoustic 

and  hydrodynamic  errors,  we  present  Fig.  3.  The  same  parameter 

values  are  used  as  in  Fig.  2,  except  e and  y are  increased  to 
1/2 

A (roughly  ten  per  cent) . For  any  uQ  and  fixed  frequency, 
the  allowable  upper  limit  for  range  is  at  least  as  large  as  the 
corresponding  value  in  Fig.  2,  as  a result  of  the  significantly 
larger  error  tolerances.  However,  the  boundary  for  consistent 
inclusion  of  hydrodynamic  effects,  represented  by  the  lower 
line,  has  moved  upward.  Thus,  the  smaller  range  and  lower 
frequency  region  is  voided,  because  resolution  of  the  environmental 
effects  is  diluted  by  larger  permissible  errors.  For  one 
comparison,  = 1 cm  sec  ^ and  f^  = 100  Hz  allows  values  of 
of  approximately  30  to  1000  Km. 
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The  cases  just  discussed  are  meant  to  be  representative 
of  typical  conditions.  Other  parameter  values  or  error  tolerances 
may  require  new  interfacing  trapezoids.  In  any  case,  the  role 
of  systematic  model  interfacing  is  essential,  as  illustrated 
here  in  constraining  acoustical  operating  conditions.  Under 
the  assumptions  of  our  environmental-acoustic  model,  the  constraints 
are  necessary  and  sufficient  for  its  application.  Other  investigators 
have  developed  different  restrictions  on  the  application  of  ray- 
theoretic  models,  based  on  different  assumptions.  For  a recent 
example.  Ref.  4 illustrates  limitations  due  to  scattering  of 
rays  from  sound-speed  fluctuations,  for  a particular  internal- 
wave  model  and  ray  type.  However,  interfacing  requirements 
supply  additional  necessary  conditions  that  must  be  satisfied  by 
the  acoustical  operating  conditions.  To  the  authors'  knowledge, 
this  interfacing  process  has  been  overlooked  previously. 

Provided  the  consistency  conditions  hold  with  y > A,  further 

simplifications  in  the  order  one  acoustic  quantities  are  possible. 

It  is  consistent  to  neglect  all  terms  of  order  y in  the  integrand 

of  the  phase  <J>,  not  only  those  discarded  in  the  derivations  of 

expressions  for  c and  u.  First,  we  note  from  Eq.  (14b)  that 

c = 1 + 0(A).  Since  we  have  restricted  attention  to  SOFAR 
s 

1/2 

rays,  we  obtain  from  Eq.  (10a)  that  both  £ and  6 are  order  A 
Therefore,  to  order  y,  r > reduces  to: 

X 

<f>  = - F f (c  + u^  + 0w)dR  , (18a) 

where  u^  is  the  component  of  the  horizontal  velocity  in  the  ng 
direction,  and  w is  the  vertical  velocity  component.  Further, 
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suppose  the  oceanic  flow  is  treated  as  incompressible.  Then  by 

scaling  the  continuity  equation,  it  follows  that  0(w/u^)  = 6 = D/L, 

where  L is  a horizontal  dynamic  length  scale.  If  for  example  L 

15  21 

is  the  Rossby  radius,  ' then  the  parameter  6 is  the  same  order 
of  magnitude  as  A.  It  follows  that  w is  consistently  neglected 
in  Eq.  (18a) : 


= - F 


A 

I 


(c  + u ) dR  . 


(18b) 


We  note  that  in  the  process  of  evaluating  or  approximating  <j), 
it  may  be  appropriate  to  neglect  terms  of  order  greater  than  the 
previously-specified  y.  In  this  event,  y should  be  replaced  by 
the  new  larger  value,  and  terms  of  this  order  should  be  eliminated 
ab  initio  in  the  derivations  of  expressions  for  u and  c. 

Similarly,  if  also  e >_  A,  it  follows  from  Eqs.  (8)  that  the 
order  one  amplitude  A,  Eq.  (lid) , is  simplified  by  neglecting 
all  terms  of  order  A in  its  evaluation.  As  a final  example, 
we  note  that  models 

for  c and  u often  depend  on  the  Brunt-Vaisala  frequency. 

Consistency  implies  that  values  for  the  square  of  the  Brunt- 
Vaisala  frequency  are  obtained  from  the  derivative  of  the  static 
sound  speed,  using  Eqs.  (14b)  - (14d) : 


N2  (z)  = v-1  [ (y-v)-s' (z) ] 


(18c) 


We  note  that  Eq.  (18c)  explicitly  demonstrates  the  connection 
between  two  essential  environmental  and  acoustical  quantities, 
the  Brunt-Vaisala  frequency  and  static  sound-speed  distribution, 
respectively. 


in  this  section,  we  apply  the  acoustic  model  of  Sec.  I 

. 2 
to  derive  expressions  for  the  variance  a of  the  phase  for  a 

SOFAR  ray  propagating  through  a stochastic  field  of  internal 

waves.  An  approximation  to  this  quantity  is  derived  and 

investigated  in  Sec.  IV.  The  ray  phase  variance  is  of 

fundamental  interest  for  at  least  two  reasons.  First,  from 

Eqs.  (11),  the  internal-wave  field  affects  only  the  phase  of  a 

ray  at  the  receiver.  Since  the  mean  of  the  phase  perturbations 

due  to  internal  waves  is  zero,  the  phase  variance  is  the  simplest 

statistic  influenced  by  the  internal  wave  field.  Second,  it  is 

necessary  and  sufficient  to  determine  the  expected  intensity 

under  the  conditions  of  Eq.  (12b). 

The  calculation  of  ray  phase  variance  has  been  treated  by 

f , 3«.Q 

other  investigators.  However,  the  development  in  this 

and  the  following  section  differs  from  their  approaches  in  at 
least  one  of  six  features:  (1)  an  "honest"  stochastic  approximation,27 
as  formulated  in  Sec.  I,  is  used  for  the  acoustic  model;  (2) 
we  are  specifically  interested  in  non-axial  SOFAR  rays,  as  opposed 
to  just  the  special  case  of  the  axial  ray;  (3)  we  develop  our 
results  for  a general  static  sound -speed  distribution;  (4)  a 
consistent  internal-wave  model  is  used,  including  characteristics 
such  as  vertical  mode  structure;  (5)  our  formulation  is  independent 
of  the  specific  internal-wave  spectral  models  used;  and  (6)  the 
effects  of  the  internal-wave  velocity  field  are  incorporated 
explicitly.  i 


l 
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As  a model  for  the  internal-wave  field,  we  use  the  formulation 
developed  in  Ref.  15.  The  static  sound  speed  is  deterministic, 
while  the  velocity  u and  order  Mach  number  sound  speed  c are 
stochastic  processes.  The  complete  specifications  of  the 
ensembles  for  these  processes  are  contained  in 
Ref.  15.  Principal  properties  include  zero 

expected  values  and  wide-sense  stationarity  in  the  variables 
x,  y,  and  t but  not  in  z.  We  wish  to  impose  here  the 
restrictions,  such  as  continuously-refracting  mean  rays,  regular 
ray  perturbations,  etc.,  under  which  the  ray  phase  and  amplitude 
formulas,  Eqs.  (11),  were  derived.  In  particular,  in  Sec.  I 
sound-speed  and  current  fluctuations  are  assumed  time  independent; 
but  this  is  not  valid  for  the  internal-wave  field.  However,  the 
dynamic  time  scale  for  internal  waves  is  about  four  hours,  and 
the  correlation  time  of  the  stochastic  internal-wave  field  is 
approximately  1.5  hours.15  Since  these  time  scales  are  much  longer 
than  acoustic  travel  times,  the  internal-wave  field  may  be  treated 
as  quasi-steady,  and  any  time  variations  in  it  may  be  studied  by 
changing  values  of  parameters  by  which  it  is  characterized. 

Also,  the  horizontal  and  vertical  correlation  lengths  of  the 
internal-wave  field  are  depth-dependent  and  are  illustrated  in 
Figs.  3 and  4 of  Ref.  15.  In  view  of  those  results,  we  continue 
to  use  the  length  scale  Xq  = D,  as  specified  in  Sec.  II. 

Under  above  assumptions,  the  ray  phase  variance  is  the 
variance  of  the  order  one  phase  of  Eq.  (lib).  However,  for  the 


i 


0 


internal-wave  formulation  of  Ref.  15,  the  parameter  6 is  of 
order  A.  Hence  from  Sec.  II  with  y >.  A,  it  follows  that  <)>  is 
given  by  the  simpler  expression  Eq.  (18b).  The  phase  variance  is 
then 


2 

o = < 


X X 


4»2  > = F2  j | cP(R1,R2)dR1dR2 


(19a) 


'0  0 

where  the  integrand  cP  is 


;£)  - < (c+u^  )|  (c+u^  > 

' (R^,£»ns) 


(R2 »£»n  ) 


(19b) 


We  note  from  Eq.  (19b)  and  the  definitions  of  R,£,n  in  Sec.  I 
that  is  the  autocorrelation  of  the  random  function  c + u 

n 

between  the  spatial  points  [x^  + R^  cos  y^  + R^  sin  T)gr 

z^R^S)]  and  [x^  + R2  cos  ng*  + R2  sin  n8,  z2(R2,£)]. 

A considerable  simplification  of  J?  is  achieved  by  assuming 

that  the  internal  wave  field  is  horizontally 

isotropic,  i.e.  that  there  exists  no  preferred  horizontal 

direction.  In  this  case,  it  follows  from  Eq.  (11)  of  Ref.  15  that 

the  cross-correlation  between  the  sound  speed  and  velocity 

fluctuations  is  zero.  Further,  since  the  mean  rays  z = z(R,£) 

are  cylindrically  symmetric,  it  can  be  shown  with  considerable 

algebra  that  is  independent  of  n . To  evaluate  , we  set 

s 

n = 0 to  obtain 
s 

Jt  R,.R2  *•£)  = < cl  cl  > + < u|  u|  > , (19c) 

'(Rj/O  1 (R2#5)  '(Rj/5)  '(R2'5) 


.are 


] 
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where  u at  n = 0 equals  u.  Thus,  the  phase  variance  splits 

n 5 

into  individual  contributions  from  autocorrelations  of  sound- 
speed  and  velocity  fluctuations.  This  was  observed  in  Ref,  1 
for  different  internal-wave  and  acoustic  models  (see,  for  example, 
features  (1),  (4),  and  (5)). 


It  is  useful  to  consider  Eq.  (19c)  in  connection  with 

reciprocal  sound  transmissions,  i.e.  two-way  propagation  between 
-►  -►2 

hydrophones  at  x^  and  x ^ . Let  the  phase  for  propagation 
in  the  direction  from  x^  to  x^  minus  that  from  x^  to  x^>  , 
along  the  same  mean  ray  at  the  same  time,  be  defined  as  the 


excess  phase  <J>  . Using  Eqs.  (11a)  - (11c),  we  obtain 
£ 


♦E  = - 2p  r *•» 

•'o 


dR  , 


(20a) 


and  note  that  <J>E  depends  on  fluctuations  only  through  the 

2 

velocity.  The  variance  of  <{>  is  denoted  by  £ and  is  given  by 


£2  = 4F: 


X X 


/f  < u|  u|  > dR. dR  . 

J (R. * C)  ' (R_»C)  1 2 

0 0 1 2 


(20b) 


Similarly,  we  define  the  additive  phase  <f>  as  the  sum  of  the 

A 

two  phases  used  to  define  <|>E.  It  follows  from  Eqs.  (11a)  - (11c) 
that 


<(>  = 2wT  - 2F  f 

Jo 


c dR  . 


(20c) 
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The  additive  phase  depends  on  fluctuations  only  through  the 

2 

sound  speed,  and  its  variance  Q is 


2 2 
CL  = 4f2 


C|  > dR.dR  . ( 2 Od ) 

a 1 (R2»C)  1 2 


From  Eqs.  (19a),  (19c),  (20b),  and  (20d)  the  total  phase 
variance  of  a ray  is  related  to  the  variance  constituents 
as  follows: 


a2  = 1/4  IQ}  +£2)  . 


(20e) 


Thus,  reciprocal  transmissions  permit  interpretation  of  the 
variance  constituents.  Furthermore,  suppose  an  acoustic 
experimental  configuration  is  available  so  that  the 
constituents  OL  and  £ are  each  individually  measurable. 
Then,  comparisons  of  their  values  may  be  used  to  estimate 
internal-wave  parameters,  as  will  be  illustrated  briefly  in 
Sec.  IV. 

The  autocorrelations  appearing  in  Eqs.  (20b)  and  (20d) 
are  obtained  from  Eqs.  (18)  of  Ref.  15: 

< c(Ri,z1)c(R2,z2)  > = (21a) 

2v2n"1N3/2(z1)N3/2(z2)  • 

OD 

H(j/j*)sin  [ jwG  < *x)  1 sin[jirG(z2)  ] • 

j-1 

B(w)  (l-u)"2)J0[fi(R2-R1)a(j,w)]dw 
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and 


< u(R1,z1)u(R2,z2)  > = 


(21b) 


n:1N1/2(z1)N1/2(z0)  • 


J) H ( j/j  *) cos  I j ttG  ( z ^ ) ) cos  [ j nG ( z2)  ] 
j=l 


(f 


'-2 


B(o>)  ^(1+u)  )Jq(6  (r2-ri) a ( j *w) ] 


A A 

-(!-«"  )J2[6  (R2~R1) a ( j »w) 


l]dw  | 


In  Eqs.  (21),  N.  = I N(z)dz  is  the  depth  average  of  the 

J0 

dimensionless  Brunt- Vaisala  frequency  N(z).  The  function 

N(s)ds  is  the  average  cumulative  Brunt-Vaisala 
frequency.  We  assume  that  the  internal -wave  areal  energy 

A 

density  spectrum  E(w,j)  is  equal  to  the  product  of  a continuous 
marginal  frequency  spectrum  B(a>)  and  a discrete  marginal  mode 

A 

number  spectrum  H ( j/  j * ) . The  internal-wave  frequencies  ui  are 
scaled  such  that  one  corresponds  to  twice  the  earth's  local 
normal  component  of  angular  velocity  (the  inertial  frequency) . 
The  vertical  mode  numbers  j of  the  internal-wave  field  are 
scaled  by  the  similarity  parameter  j*,  an  effective  mode-number 
bandwidth.  The  marginal  spectra  are  normalized  as  follows: 

00  00 

^ B (o>)  dio  = £ H ( j/  j*)  = 1 , (21c) 

1 j=l 

so  that  the  total  dimensionless  expected  energy  per  unit  area 


is  one.  This  condition  is  dimensionally  satisfied  by  choosing 
the  characteristic  velocity  UQ  * (Ep/ppD)*/2,  where  Ep  is 
the  total  dimensional  expected  energy  per  unit  area.  The 
value  Up  = 2.72  cm  sec-*  corresponds  to  eq  = 3.82  x 103  Jm-2. 
The  wave  numbers  a of  the  internal  wave  field,  appearing  in 
the  arguments  of  the  order  n Bessel  functions  J^,  are  related 

A 

to  a)  and  j via  the  dispersion  relation 

-1  ~ 2 1/2 

a - jirNp  (u>  -1)  ' . (21d) 

The  dimensional  wave  numbers  are  a/L,  where  L is  the  Rossby 
radius,  so  that  the  parameter  6 multiplies  (R2~R^)a  after  the 
range  is  scaled  by  D. 

Upon  substituting  Eqs.  (21)  into  Eqs.  (20b)  and  (20d) , we 

A O 

obtain  lengthy  expressions  for  the  variances  CL  and  £ . in 

general,  simpler  analytical  formulas  are  not  available  for 

these  integrals,  and  numerical  or  asymptotic  methods  are 

employed  for  their  evaluation.  There  are  two  classes  of  SOFAR 

rays  for  which  the  integrals  can  be  analytically  approximated. 

Let  the  parameter  tc  = x Max  |d*/dR|  , the  ratio  of  transmission 

R 

range  to  the  scale  length  of  the  mean  ray.  From  Eq.  (10b) , 
k = x tan  IO  • Where  6 is  the  ray  angle  at  the  SOFAR  axis. 
The  first  class  of  rays  is  distinguished  by  k<<  1.  included  in 
this  case  are  the  axial  ray  (k  * 0)  and  "near  axial"  rays 
with  relatively  small  depth  variation  with  range  or  values  of  x 
near  one.  The  range  integrals  for  the  axial  ray  z * a,  where  a 
is  the  SOFAR  axis  depth,  have  been  evaluated  previously.6'7 
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A second  class  of  rays  consists  of  those  rays  for  which  k >_  1, 
and  is  discussed  in  Sec.  IV.  Prom  Eq.  (10a)  this  condition  is 
equivalent  to 

X > [ca  (cos5)/C>|  ]-{l  - lca(cosO/c^  iV172  (22a) 

where  c = c (a) . it  follows  that  a sufficient  condition  for 
a s 

1/2 

the  validity  of  Eq.  (22a)  is, to  order  A ' , 

x > [2  (cj  -ca)]"1/2  . (22b) 

He  note  that  Eq.  (22b)  is  typically  more  restrictive  than  the 
condition  x 1 1 assumed  in  Sec.  n.  For  a numerical  example, 
if  the  dimensional  source  and  axial  sound  speeds  differ  by 
30  m sec  * , and  if  the  transmission  ranges  are  greater  than 
about  25  km,  all  SOFAR  rays  satisfy  the  condition  < >_  1. 

IV.  PHASE  VARIANCE  FOR  tc  > 1 

In  this  section,  we  approximate  asymptotically  the  expressions 

for  the  variance  constituents  CL  and  <£,  . The  phase  of  the 

sinusoidal  terms  appearing  in  Eqs.  (21a)  and  (21b),  and 

consequently  in  the  integrands  of  the  variance  constituents, 

-1  fz 

are  of  the  form  jir  Nn  1 N(s)ds.  In  the  derivations  of  Eqs.  (21), 

0 Jo 

it  has  been  assumed^  that  ju/Ng  is  a large  parameter.  Therefore, 

28 

the  method  of  stationary  phase  may  be  used  to  simplify  the 
integrals  over  range  (R^  and  r2)  of  the  variance  constituents. 


The  stationary  points  are  values  of  and  R2  for  which  dz/dR 
vanishes,  i.e.  turning  points  of  the  SOFAR  ray.  It  can  be 
shown  that  this  method  is  applicable  if  x (dz/dR)  is  at  least 
order  one,  with  respect  to  the  large  parameter  jTr/NQ»  at  some 
values  of  r on  the  ray.  Thus,  it  is  sufficient  to  require 
< ^ 1.  Henceforth  in  this  section,  we  assume 
the  more  restrictive  condition  Eq.  (22b). 

We  emphasize  that  it  is  the  inclusion  of  the  j-dependent 
modal  structure,  c.f.  feature  (4)  in  Sec.  Ill,  which  makes 
application  of  stationary  phase  feasible.  Other  recent 
4 5 

investigations  ' have  used  a type  of  turning  point  method,  an 
apex  approximation,  but  without  apparent  asymptotic  justification. 
Thus,  the  consistent  formulation  of  our  internal-wave  model  provides 
a theoretical  justification  for  an  apex  approximation.  Moreover, 
different  numerical  predictions  are  obtained,  as  will  be 
illustrated  later.  In  addition,  use  of  stationary  phase  extends 
the  validity  of  the  approximation  to  include  turning  points 
below  the  SOFAR  axis,  not  just  to  those  above  the  axis.4 

Applying  stationary  phase  to  first  the  Rj^  and  then  the  R2 
integral  in  the  expressions  for  CL  and  , we  obtain 


CL2  = 


(8v2F2/A) 


K 

£ 

i=i 


K . 

£ 

m=l 


7!  [H(  j/j*)/jj  sin  [ jirG(z  ) + ir/4]  sin  [jirG(z  ) + it/4]  • 

x,  — m — 

j=l 


{(  B(w)  (l-i‘2)J  [6(R  -R.) 


a(j,w)]dw}l 


(23a) 
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K K 


& = <4F2/4>  £ £ 

£=1  m=l 


Is'  u£)s'  («B)  r1/2  • 


Y cos  [jirGtZj, ) + ir/4]  cos  [ jirG  (zm)  + ir/4] 

j=l 


B (01)  { 


B (oj)  (l+u"2)J0[6(Rm-RJl)  a ( j ,u) ] 


- (l-w_2)J2[6(Rm-Ra) 


a(j,u>) 


doj 


(23b) 


In  the  derivation  we  have  used  the  expression  for  the  ray 
2 2 

curvature  d z/dR  = - As1 (z) , which  follows  from  Eqs.  (10b) 
and  (14b)  and  consistent  neglect  of  0(A)  terms.  In  Eqs.  (23), 
the  integer  K is  the  total  number  of  ray  turning  points  for 
the  range  x*  The  ith  turning  point  is  located  at  range  R^  and 
depth  zA,  labeled  so  that  R^<R^<  • • • < R^ . The  choices  of 
signs  in  the  sinusoidal  arguments,  jirGfz^  + ir/4,  are  determined 
as  positive  if  s'  (z^  < 0 and  negative  if  s'  (z^  > 0. 

Suppose  a ray  is  regarded  as  divided  into  segments;  the 
first  consists  of  the  portion  between  the  source  and  the  first 
SOFAR-axis  crossing,  the  second  between  the  next 
consecutive  axis  crossings,  etc.,  and  the  last 
from  the  final  axis  crossing  to  the  receiver.  From 
stationary  phase,  the  main  contribution  to  CL*  (or  <£/ ) from 
each  ray  segment  comes  from  a neighborhood  of  its  turning  point. 
If  an  initial  or  terminal  segment  has  no  turning  point,  its 
contribution  is  absent.  Therefore,  in  view  of  Eq.  (20d)  (or 
Eq.  (20b)),  we  interpret  the  Jt  / m terms  in  the  sums  of  Eq. 

(or  Eq.  (23b))  as  contributions  from  the  autocorrelations  of 


(23a) 
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sound-speed  (or  velocity)  fluctuations  between  neighborhoods  of 
turning  points  in  distinct  ray  segments.  For  physically 
interesting  models  of  B and  H (for  example.  Refs,  lo  and  29), 
the  autocorrelations  of  sound  speed  and  velocity  decrease 
significantly  in  magnitude  for  separations  in  depth  of  about 
2 km  and  in  range  of  about  15  km.15  Thus,  provided  the 
respective  horizontal  and  vertical  distances  between  turning 
points  are  greater  than  these  correlation  lengths,  the  £ ^ m 
terms  in  Eqs.  (23)  are  negligible  compared  to  the  £ = m terms. 
This  assumption  is  typically  satisfied  for  SOFAR  propagation  in 
the  deep  ocean. 

In  addition  to  neglecting  £ ^ m terms  in  Eqs.  (23) , we 
assume  that  each  mean  SOFAR  ray  has  exactly  two  distinct  turning 
depths.  He  denote  the  turning  depth  above  (or  below)  the  axis 
by  z+(or  z ),  where  s' (z+)  < 0 (or  s' (z  ) > 0),  and  we  note  that 
s(z+)  « s ( z ).  Also,  we  let  k+  (or  K~)  be  the  number  of  turning 
points  above  (or  below)  the  axis.  It  is  convenient  to  note  that 

s' (a)  * 0 * (y-v)  - vN2(a)  , (24a) 

which  follows  from  Eq.  (18c)  and  the  definition  of  the  SOFAR  axis 

2 + 2 

z = a.  Assuming  the  typical  situation  of  N (z  ) > n (a)  and 
2 - 2 

N (z  ) < N (a),  it  follows  from  Eqs.  (18c)  and  (24a)  that 
i 2 - 7 

| s ' ( z )|  = +v[N  (z  ) - N (a)]  . (24b) 

2 2 

We  also  define  the  ratio  n (z)/n  (a) 

+ 

8(z  ) £ !• 


= p(z),  so  that  typically 


r 
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Using  Eq.  (24b),  we  obtain  compact  expressions  for  Eqs.  (23), 
2 


A.  - d q ik  qc  + k-q;i 


(25a) 


and 


£ 2 - <£  * ikX  + KX]  > 

where  we  have  introduced  the  following  notation: 
2 


(25b) 


&Q  = 4vF2(1-I)/ (AJa)  , 

(26a) 

£0  = 2F2(1+I)/[(v-v)AJ*]  , 

(26b) 

Q*  * 8(8-1)  "V  +J[G])  1 , 

c 1 z«=z 

(27a) 

Q~  - 6(1-6)  IGJ  ) | _ , 

z— z 

(27b) 

Q*  = (8-l)_1 IGJ ) 1 , 

1 z*z 

(27c) 

q"  = (1-8)  _1(1+J  [G] ) 1 

' Z=Z 

(27d) 

I - f [B(w)/u2]dw  , 

(28a) 

*1 

OO 

J*1  - £ H ( j/j*)/j  , 

(28b) 

j=l 


and 


oo 

J>lG(z)]  = J*  £ [H(j/j#)/j]  sin  [2ir jG(z)  ] . 

j=l 


(28c) 


The  variance  constituents  Eqs.  (25)  are  each  proportional  to 
normalization  factors,  Eqs.  (26),  which  in  turn  depend  on  two 


f\ 
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characteristic  parameters  of  the  internal-wave  field,  Eqs.  (28a) 

and  (28b) . The  number  of  upper  and  lower  turning  points  appear 

in  Eqs.  (25)  multiplied  by  Q factors,  Eqs.  (27)  . We  note  that 

typically  SOFAR  rays  for  relatively  large  source-receiver  separations 

have  more  turning  points  than  those  for  relatively  small  ranges. 

Thus,  the  linear  dependence  of  and  (£  2 on  K+  and  K~  reflects 

a general  tendency  for  their  increasing  values  with  range,  but 

their  exact  variation  with  x is  a complicated  function  of  ray 

geometry.  The  quantity  Q+  represents  the  contribution  due  to 

c 

sound-speed  fluctuations  at  one  upper  turning  point  z = z+,  and 
the  other  Q's  have  analogous  interpretations.  These  quantities 
in  turn  depend  on  B and  the  function^,  Eq.  (28c). 

The  only  dependence  on  B(iii)  is  through  the  parameter  I, 

Eq.  (28a).  in  Ref.  15,  it  is  shown  that  the  fractional  amounts 
of  expected  potential  and  kinetic  energies  per  unit  horizontal 
area  of  the  internal-wave  field  are  equal  to  (l-l)/2  and  (l+i)/2, 

respectively.  Thus,  the  relative  contributions  of  the  spectrum 

^ 2 
B (to)  to  the  sound-speed  dependent  6C  and  velocity  dependent 

are  given  by  the  fractions  of  expected  potential  and  kinetic 

A 

energies,  respectively.  Since  0 <_  I <_  1,  the  influence  of  B(w) 

. , . 2 2 

is  to  increase  the  contribution  to  o of  , relative  to  that  of 
^2.  For  the  choice  B(<i>)  = 2 (iru)  ^(oj2-1)  the  ratjo  of 

expected  kinetic  to  expected  potential  energy  is  three  to  one. 

The  dependence  of  and  £,2  on  H(j/j*)  enters  via  the 
combinations  JA  and  d , defined  in  Eqs.  (28b)  and  (28c), 
respectively.  We  next  consider  each  of  these  in  turn.  We  note 
that  interpreting  H(j/j#)  as  a discrete  probability  density  for 
the  vertical  mode-number  energy  distribution  leads  to 
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J4  = <l/j  >.  As  an  example,  for  Cairns'  model  of  H(j/jA), 


h ( j / j * ) = (j2  + j2)  V£  (j2+j2)-1  » 

j=l 


(29a) 


we  obtain 


- £ h +j*i  / 2 ho2^)]-1  . 
j-i  j-i 


(29b) 


In  the  limit  -►  0,  Eq.  (29b)  reduces  to  = c(2)/c(3)  ~ 1.37, 
where  c (n)  is  the  Riemann  zeta  function.  For  large  values 

of  j , it  follows  that30 

* 


J,  ~ (irj*  - 1)  / (2y  + In  j3)  , 


(29c) 


where  y * 0.557  is  the  Euler  constant,  and  the  approximation 
from  Eq.  (29c)  is  excellent  for  j#  4.  The  value  of 

JA  increases  with  j#,  and  is  less  than  j * for  j > 2.  From  the 
appearance  of  the  factor  in  Eqs.  (26),$ 3 and  decrease 

as  j,  increases.  This  behavior  is  due  to  increased  mode 
interaction  and  "cancellation" . 

The  vertical-structure  function  j^|  (G) , Eq.  (28c) , is  analogous 
to  the  function  S(G), 


S<G)  = ^ H(j/j*)  cos  (2 jirG]  , 

j=l 


(30a) 


which  arises  in  the  internal  wave  model  to  apportion  relative 
amounts  of  expected  potential  and  kinetic  energies  throughout 


..  — ..  .".T  • 


I 


0 

n 
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the  depth. 


15 


i 


4 


The  relation  between  & and  S is 
G 


(G)  = 2ttJ 


* J S(g)dg 
0 


(30b) 


Whereas  the  spectral  component  B(co),  through  the  parameter  I, 
reflects  the  relative  expected  energy  contributions  to  variance 
constituents  in  a depth-independent  manner,  H modifies  those 
contributions  through  the  depth-dependent  behavior  of  ^ . This 
result  explicitly  shows  an  effect  of  our  inclusion  of  internal- 
wave  vertical  structure,  a consequence  of  feature  (4)  of  Sec.  ill. 
From  Eqs.  (28b)  and  (28c) , we  note  that  for  any  model  of 

H(j/j*)  > 0:  (a)  | ^(G)  | <_  1 ; (b)  ^(G)  = 0 for  G = 0,  1/2,  1; 

(c)  y4(G)  is  odd  about  G = 1/2;  and  (d)  ^(G)  has  a maximum 

between  G = 0 and  G = 1/2,  equal  to  the  negative  of  its  minimum, 

which  is  between  G = 1/2  and  1. 

For  Cairns'  model  of  H(j/j  ),  Eq.  (29a),  > ^ reduces  to 


i 


(G)  = ttJ* 


(1-2G)  sinh(  j*ir)  -sinh  [ j*ir  (1-2G)  ] 
j^TT  cosh(j^Tr)  - sinh(j^TT) 


(31a) 


In  the  limit  j -*■  0,  we  obtain 


^l(G)  = t(l-2G)  (G-G2)] 


(31b) 


In  Fig.  4 , we  plot  ^)[G(z)]  from  Eqs.  (31)  versus  depth,  using 
the  exponential  Brunt-Vaisala  frequency  model  N(z)  = 3/2  exp(-4z) 
and  the  values  j#  = 0,  3,  and  7.  Three  has  been  suggested  for 
j#  in  Ref.  29  , but  other  values  (see,  for  example,  Ref.  5) 


15 
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may  be  appropriate,  depending  on  the  oceanic  region.  For 
reference,  both  the  SOFAR  axis  depth,  a = .26  from  Eq.  (24a), 
and  the  depth  b = .17  for  which  G(b)  = 1/2  are  shown  dashed. 

The  values  y = 11.6  and  v = 9.1  in  this  example  are  obtained 
from  Eqs.  (14b)  - (14d)  by  requiring  that  N(z)  produce  surface 
and  axis  sound  speeds  equal  to  1540  and  1490  m sec-1,  respectively 
in  an  ocean  with  D = 5 km  and  equal  bottom  and  surface  sound 
speeds.  We  note  that  the  amplitude  3/2  of  the  exponential  of 
N(z)  has  no  influence  on  either  b or  p^[G(z)],  but  it  is  necessary 
to  determine  the  SOFAR  axis.  The  curves  in  Fig.  4 possess  the 
general  characteristics  of  the  previous  paragraph.  Their 
asymmetry  in  z results  from  the  depth  dependence  of  G(z).  The 
situation  a > b shown  in  Fig.  4 is  in  fact  typical.  Thus, 
typically  there  exist  two  types  of  upper  turning  points, 
corresponding  to  0 < z+  < b for  which  ji  >_  0,  and  b < z+  < a, 
for  which  0.  In  addition,  from  Fig.  4 , we  observe  that  as 

j*  increases,  the  maximum  value  of  |^[G(z)]|  decreases.  This 
result  is  in  concert  with  the  cancellation  effect  discussed  in 
connection  with  Eq.  (29b) . 

We  next  investigate  the  behavior  of  the  turning-point 
contributions,  the  Q factors  of  Eqs.  (27).  We  note  that  they 
depend  only  on  the  specification  of  H(j/j*),  N(z),  a,  and  z^  , 
through  the  functions  g and  d • We  plot  Eqs.  (27)  versus 
dimensionless  upper  turning-point  depth  z+  in  Fig.  5 , for  j#  * 3 
and  the  same  example  functions  and  other  parameter  values  used 
in  Fig.  4.  Functions  <£  (or  Q^)  are  shown  solid  (or  dashed). 

The  Q functions  are  evaluated  at  the  lower  turning-point  z- 
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corresponding  to  any  z (such  that  s(z  ) = s(z  ))  for  the  given 

mean  sound-speed  profile.  The  vertical-structure  function 

affects  the  Q's  at  all  values  of  z+.  its  most  important 

influence  is  on  Q+  and  Q for  z+  between  0 and  0.1,  at  which 

c u 

depths  enhances  Q*  and  detracts  from  Q+.  The  sharp  change 

in  value  of  Q+  near  the  surface  is  a particularly  obvious 

feature  due  to  d , c.f . Pig.  4 . The  disparity  in  values  of  Q+ 

and  Qu  near  z+  = 0 results  from  the  larger  mean  sound-speed 

gradient  there  as  compared  to  that  near  z~  = 1.  The  small  value 

of  Q-  near  the  surface  is  due  to  the  nearly  negligible  value  of  6 
c 

near  the  bottom.  Finally,  as  z+  and  z approach  z = a,  the  sound- 

speed  gradient  and  ray  curvature  approach  zero,  and  consequently 

0 approaches  one,  causing  the  increase  in  all  four  curves 

shown  near  the  bottom  of  Fig.  5 . However,  our  approximate 

formulas  are  not  valid  for  z+  too  close  to  the  axis,  because  of 

our  use  of  the  stationary-phase  method.  In  particular,  given  a 

nondimensional  range  Eq.  (22b)  provides  a condition 

_2 

(c^  - ca)  >_  1/2  x to  be  satisfied  by  the  minimum  source 
sound  speed.  This  in  turn  is  a condition  on  the  maximum  source 
depth,  which  is  an  upper  bound  on  z+ . 

From  Fig.  5 , we  observe  that  the  Q's  are  order  one  in 
magnitude.  Thus,  the  normalizations  of  the  variance  constituents, 
Eqs.  (26),  provide  appropriate  scales  for  individual  turning- 
point  contributions.  Suppose  now  that  the  brackets  in  Eqs.  (25) 
are  comparable.  The  relative  significance  of  sound-speed  and 
current  effects  to  the  phase  variance  is  determined  by  the 
relative  magnitudes  of^.*  and  . Let  the  ratio  r be  defined  by 


D 
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r = « q/Oo  = 2vtvv)  (1-D/U+I)  • 


(32) 


Then,  T>>1  (or  r <<  1)  corresponds  to  sound-speed  fluctuations 

dominating  current  fluctuations  (or  vice  versa) . The  value  of 

T varies  with  v,  y,  and  I.  With  respect  to  the  I variation, 

we  note  the  earlier  discussion  that  r is  proportional  to  the 

ratio  of  expected  potential  and  kinetic  energies.  The  limiting 

case  I = 0 (or  1)  represents  internal-wave  energy  concentrated 

at  the  inertial  (or  far  from  the  inertial)  frequency.  For  1*0 

there  are  equal  amounts  of  expected  energies,  while  for  1*1  only 

kinetic  energy  is  present.  Thus,  for  the  v and  y numbers  used  in 

Figs.  4 and  5 and  values  of  i * 0,  0.5,  and  1,  Eq.  (32)  yields 

respective  values  of  r = 45.5,  15.2,  and  0.  We  conclude  that  the 

relative  significance  of  sound-speed  and  current  effects  is 

strongly  dependent  on  a parameter  characterizing  the  frequency 

distribution  of  internal-wave  energy.  The  above  results  are  a 

consequence  of  our  general  internal-wave  model  formulation,  c.f. 

features  (3)  through  (6)  of  sec.  Ill,  and  have  been  overlooked 

1,2 

xn  previous  considerations  of  internal-wave  current  effects. 

There  are  obvious  exceptions  to  the  usefulness  of  Eq.  (32) 

to  determine  the  relative  significance  of  sound-speed  and  currents 

on  phase  variance.  Suppose  a ray  has  a single  turning  point 

located  deep  below  the  axis.  Then,  current  fluctuations  are 

inevitably  dominant,  since  Q~  is  nearly  zero  as  in  Fig.  5 . For 

c 

reciprocal  transmissions,  the  additive  phase  variance  is  always 
sound-speed  dominant,  and  the  excess  phase  variance  always 
current  dominant.  As  an  example  of  exploiting  reciprocal 
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transmissions,  consider  a near-grazing  ray  with  |K~  - K+|/  K+  <<  1. 

It  follows  from  Eqs.  (25),  (27),  and  (28)  that  $2/  £2  is 

approximately  r.  Since  and  £,2  are  individually  observable, 

the  value  of  I could  be  estimated  if  values  for  y and  v are  known. 

To  close  this  section,  we  compare  our  phase-variance 

predictions  with  the  only  other  calculation  for  a non-axial  ray 

known  to  us  in  the  literature,  i.e.  Eq.  (7.4.8)  of  Ref.  4.  That 

formula  corresponds  to  the  particular  case  in  our  results  of 

ignoring  (£,2  in  Eq.  (20e)  and  Q in  Eq.  (25a) . It  can  be  shown 

c 

that  using  the  parameter  values  of  Ref.  4,  our  phase-variance 
prediction  differs  by  the  factor  1 + >4tG(z+)].  From  Fig.  5,  it 
follows  that  for  turning  points  between  the  surface  and  z «=  b, 
i.e.  z * = 850  m,  our  prediction  is  up  to  60%  greater.  For  a 
turning  point  below  z * b,  our  prediction  is  up  to  308  less.  For 
other  parameter  values,  the  predictions  can  differ  by  either 
greater  or  smaller  amounts.  In  any  case,  this  comparison  serves 
to  illustrate  a representative  numerical  difference  arising  from 
the  features  of  our  environmental-acoustic  model. 

V.  EXPECTED  INTENSITY 

We  investigate  the  expected  value  (i.e.  the  ensemble  average) 
of  the  acoustic  intensity  in  the  presence  of  a stochastic  internal 
wave  field,  in  the  internal  wave  model,  velocity  and  sound- 
speed  fluctuations  are  random  because  of  stochastic  superposition 
coefficients.  ^ We  assume  here,  as  elsewhere  (see,  for  example. 


Ref.  13)  that  the  coefficients  have  zero-mean  normal  distributions. 
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It  follows  from  Eq.  (lib)  and  Refs.  15  and  31  that  the  order  one 

phase  of  any  ray  at  the  receiver  is  also  normally  distributed. 

Assuming  further  that  phases  of  distinct  rays  at  the  receiver  are 

mutually  independent/  the  expected  intensity,  proportional  to 
-2 

<ct  > , is  given  by  Eq.  (12b)  . Evaluation  of  this  quantity  requires 

the  amplitude,  mean  travel  time,  and  phase  variance  for  each  ray 

at  the  receiver,  along  with  the  source  frequency. 

We  specify  the  static  sound  speed  with  a bilinear  profile, 

which  has  been  employed  previously  (see,  for  example.  Refs.  32  and 

33  ) as  a convenient  mid-ocean  model.  We  take  the  surface  and 

bottom  sound  speeds  equal,  eliminating  reflected/refracted  rays. 

Furthermore,  we  assume  that  due  to  bottom  and  scattering  losses, 

surface-reflected/bottom-reflected  rays  are  of  negligible  importance 

to  the  total  received  field.  Thus,  only  SOFAR  rays  are  considered, 

and  the  ray  formulas  of  Sec.  I apply.  In  particular,  Eqs.  (10), 

(11a) , and  (lid)  provide  the  complete  specification  of  the  mean 

rays.  A further  simplifying  assumption  is  that  the  source  is 

placed  above  the  SOFAR  axis  and  the  receiver  below,  such  that 

they  are  at  equal  static  sound  speeds.  The  analysis  of  the  ray 

formulas  for  these  static  conditions  is  comprehensively  detailed 

in  Sec.  I of  Ref.  32.  We  select  the  numerical  values  of  the 

34 

bilinear  parameters  from  a recent  fit  to  deep-ocean  sound-speed 

observations:  c * 1543  m sec  \ c c = 1484  m sec  \ d = 4900  m, 

0 a o 

and  aD  = 1100  m. 

For  the  phase  variance,  we  wish  to  apply  the  formulation  of 
Secs.  Ill  and  IV.  The  principal  restriction  on  those  results  is 
Eq.  (22b),  derived  from  the  condition  k 1.  For  the  bilinear 
profile  and  source  at  nondimensional  depth  d < a,  Eq.  (22b) 
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becomes 

x>  (a/[2(l-ca) (a-d)] }1/2  . (33) 

Using  our  profile  parameters  and  dimensional  source  depths  dD  of 

100,  400,  700,  and  1000  m,  Eq.  (33)  is  satisfied  for  ranges 

greater  than  about  19,  22,  29,  and  60  km,  respectively.  As  the 

source  approaches  the  axis,  larger  ranges  are  needed  to  eliminate 

"near-axial"  rays.  However,  as  suggested  by  the  numerical 

examples  for  shallow  or  moderate  source  depths,  the  restriction 

is  not  severe  for  the  bilinear  profile.  This  results  mathematically 

-1/2 

because  Eq.  (33)  behaves  like  (a-d)  as  d -*■  a,  whereas  for 
differentiable  profiles,  the  analogous  expression  grows  faster, 
like  (a-d)-1. 

To  specify  the  variance  constituents,  Eqs.  (25),  we  first 
require  the  distribution  of  Brunt-Vaisala  frequency  n(z).  For 
our  consistent  environmental-acoustic  model,  specification  of  the 
static  sound  speed  determines  N(z),  as  discussed  in  Sec.  II.  By 
Eqs.  (14b)  - (14d) , the  bilinear  profile  leads  to  an  N(z)  which 
is  bi-constant  (two  discontinuous  vertical  line  segments) . We 
choose  the  previous  numbers  y = 11.6  and  v = 9.1,  to  facilitate 

comparisons  with  earlier  results;  also,  these  values  correspond 

* 

to  a smooth  static  profile  in  Sec.  IV  that  approximates  the 
observations  to  which  the  bilinear  profile  was  fit  in  Ref.  34. 

Thus,  we  obtain  Brunt-Vaisala  frequencies  of  4.0  eye  h”*  (or 

0.4  eye  h 1)  between  surface  and  axis  (or  axis  and  bottom).  Also, 

we  choose  the  velocity  scale  Ug  * 2.72  cm  sec-1  as  in  Sec.  III. 

The  turning-point  contributions  to  the  phase  variance  are 
specified  by  Eqs.  (27)  if  we  select  a model  for  H(j/j*),  for 
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which  we  use  Eq.  (29a).  In  Fig.  6,  we  graph  the  turning-point 
contributions  versus  dimensionless  upper  turning  depth  for  the 
bi-constant  N(z)  and  the  value  j*  = 3.  It  is  important  to  note 
that  the  curves  in  Fig.  6 enjoy  the  same  features  as  those  in 
Fig.  5 for  the  exponential  N(z),  with  one  exception.  The  curves 
in  Fig.  6 do  not  show  the  sharp  increase  as  z+  approaches  a. 

This  is  because  the  increase  results  from  profile  curvature, 
which  is  absent  from  the  bilinear  model.  Thus,  in  order  that  the 
bi-constant  and  exponential  cases  have  similar  turning-point 
contributions,  we  want  ray  upper  turning  depths  to  avoid  depths 
influenced  by  profile  curvature.  One  way  to  assure  this  is  to 
place  the  source  sufficiently  far  from  the  axis,  for  example, 
at  a dimensional  depth  of  400  m shown  on  Fig.  6 as  a horizontal 
line  (this  particular  value  was  also  used  in  Ref.  34).  All 
upper  turning  points  necessarily  lie  above  the  source  depth. 
Moreover,  comparison  of  corresponding  curves  in  Fig.  5 and  6 show 
that  the  numerical  values  above  this  source  depth  are  quite 
similar. 

Before  computation  of  the  expected  intensity,  it  is 
appropriate  to  check  an  interfacing  trapezoid  for  our  application. 

A proper  choice  for  the  error  parameter  y of  Sec.  II  is  NQ/ir=  0.10, 
which  is  the  order  to  which  the  WKB(J)  solutions  for  the  internal- 
wave  model  are  obtained  in  Ref.  15.  It  happens  that  the  WKB(J) 
solutions  are  exact  for  the  bi-constant  N(z).  However,  in  the 
stationary-phase  approximation  of  4 in  Sec.  IV,  terms  of  this 
> order  are  ignored.  Thus,  we  have  an  example  of  the  situation 
noted  just  below  Eq.  (18b),  and  our  choice  y = NQ/ir  > A is 
justified.  We  also  select  c = y for  our  tolerable  acoustic  error 
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(approximately  ten  per  cent).  From  Eqs.  (16)  and  (17),  we  may 

construct  an  interfacing  trapezoid.  We  do  not  present  its  graph, 

because  for  the  parameter  values  specified  in  this  section,  its 

three  nonvertical  sides  are  quantitatively  very  close  to  those  of 

the  UQ  = 3 cm  sec-1  trapezoid  in  Fig.  3.  Its  leftmost  (vertical) 

boundary  is  the  line  x = 4.4,  corresponding  to  equality  in  Eq.  (33) 

and  representing  a stronger  condition  than  was  required  in 

Sec.  II.  Therefore,  our  selection  of  acoustic  operating  conditions 

must  lie  within  this  trapezoid,  in  addition,  we  desire  that  the 

range  interval  lie  in  a "stability  region"  for  the  static  bilinear 
34 

profile.  Such  requirements  are  satisfied  by  the  choices 
fQ  = 100  Hz  and  RQ  between  200  and  222  km,  which  we  will  use 
for  subsequent  calculations.  We  note  that  in  this  range  interval,  four 
SOFAR  rays  intercept  the  receiver,  each  with  seven  axis  crossings. 

We  turn  to  investigation  of  the  intensity.  In  order  to 
illustrate  the  nature  of  internal-wave  effects,  we  introduce  the 
following  notation: 


=2  - £ ' 
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£A<n)A(„,)  008  1“(T(n)-T(m)))  ' 

(34b) 

n 

m 

n?*m 

and 


■ »2  + £ 

^ A(n)A(m) 

(34c) 

n 

m 

n^m 

— 2 -2 

As  mentioned  in  Sec.  I,  a (or  a ) is  proportional  to  the 

X c 

incoherent  (or  coherent)  intensity,  and  is  the  limit  of  <S2>  as 
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the  phase  variance  on  each  ray  approaches  infinity  (or  zero) . 

-2 

From  Eq.  (34b),  it  follows  that  aM  in  Eq.  (34c)  is  proportional 
to  the  maximum  value  of  the  coherent  intensity.  We  note  that 
a~2  is  the  square  of  the  sum  of  the  static  ray  amplitudes. 

Figure  7 is  a graph  versus  dimensional  range  of  Eqs.  (34) 
in  dB  relative  to  a sphere  at  dimensionless  radius  r,  as 
discussed  prior  to  Eq.  (7a).  In  addition,  <a  > from  Eq.  (12b) 
is  shown  in  dB  for  internal-wave  parameters  I = 0.5  and 

= 3.  Calculations  were  made  using  a range  increment  of  0.2  km. 
The  comparative  magnitudes  of  variations  are  as  follows.  The 
drop  in  incoherent  intensity  over  about  20  km  of  range  is 
approximately  7 dB.  The  difference  between  incoherent  intensity 
and  the  maximum  coherent  intensity  is  about  6 dB  at  each  Rq. 

The  difference  between  the  coherent  and  incoherent  intensities 
varies  from  about  -10  dB  to  6 dB,  except  for  one  deep  fade  of 
at  least  18  dB  near  212.6  km.  The  difference  between  the  expected 

and  incoherent  intensities  varies  from  about  -3  dB  to  4 dB.  We 

\ 

conclude  that  the  deviations  from  both  coherent  and  incoherent 
intensities  produced  by  internal  waves  can  be  very  significant. 
Although  at  isolated  ranges  all  three  intensities  agree,  neither 
coherent  nor  incoherent  intensities  provides  an  accurate  picture 
of  range  variations  in  the  presence  of  internal  waves.  For  example, 
the  depth  of  fades  is  curtailed  by  the  internal  waves.  We  note 
that  the  value  of  the  expected  intensity  at  each  range  lies  between 
the  coherent  and  incoherent  intensities.  From  Eq.  (12b),  this 
occurs  because  all  rays  have  approximately  the 
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same  phase  variance  for  I = 0.5  and  j = 3. 

We  now  consider  variations  in  the  expected  intensity  resulting 
from  changes  in  the  internal-wave  parameters.  The  internal-wave 

A 

energy-density  spectral  functions  B(to)  and  H(j/j*)  are  useful 

representations  of  observations.  However,  at  present  there  is 

no  certainty,  based  on  dynamical  principles,  of  the  correct  form 

or  parametric  values  for  the  spectra,^  in  any  given  ocean  region. 

We  regard  variations  in  the  parameters  j#  and  I as  modeling  certain 

types  of  changes  in  the  internal-wave  field  between  different 

oceanic  regions,  or  between  different  times  in  the  same  region. 

Indeed,  changes  in  these  characteristics  have  evolved  with  the 

10  29 

acquisition  of  new  data.  ' 

In  Fig.  8 , we  graph  the  expected  intensity  relative  to  the 

incoherent  intensity  in  dB  versus  range  for  I = 0.5  and 

j*  * 0,3,  and  7.  Thus,  the  solid  curve  in  Fig. 8 represents  the 

difference  between  the  two  solid  curves  in  Fig. 7 . We  note  the 

increased  variation  of  the  expected  intensity  about  the  incoherent 

as  j#  increases.  For  jA  = 0,  3,  and  7,  intensity  varies  by 

approximately  2.5,  7,  and  9 dB, respectively , as  RQ  changes  over  the 

interval  shown.  These  results  are  consistent  with  increased  mode 

interaction,  mentioned  in  connection  with  J*  in  Sec.  IV. 

We  show  similar  curves  with  = 3 and  I = 0,  0,5,  and  1 in 

Fig.  9 , where  the  solid  curve  has  the  same  meaning  as  in  Fig.  8. 

Again,  increased  fluctuation  of  the  expected  intensity  about  the 

incoherent  occurs  as  the  internal-wave  parameter  increases.  For 

1=0,  0.5,  and  1,  intensity  varies  by  about  3,  6.5,  and  21  dB, 

respectively,  as  changes.  The  behavior  of  the  expected  intensity 

with  I is  interpreted  by  considering  the  maximum  of  £ ^ and  as 

2 

an  appropriate  scale  for  o . From  Eqs.  (26),  it  can  be  shown 
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that  for  2v(y-v)  > 1,  this  maximum  decreases  as  I increases. 

The  fluctuations  become  very  dramatic  as  I approaches  one, 
reaching  nearly  the  size  of  those  of  the  coherent  intensity  about 
the  incoherent.  Thus,  mean  intensity  is  even  more  sensitive  to 
variations  in  I than  j*.  We  note  that  variations  with  range  in 
received  intensity  have  been  studied  previously;  recent  examples 
include  Ref.  33  and  35.  However,  to  the  authors'  knowledge,  the 
explicit  effects  of  internal  waves,  and  of  course  variations  in 
the  internal-wave  field,  on  multipath  expected  intensity  have  not 
been  considered  previously. 

VI.  SUMMARY 

The  two  principal  purposes  of  this  paper  are  to  formulate  a 
consistent  environmental-acoustic  model  for  a deep  moving  ocean  and 
to  apply  the  model  to  investigate  effects  of  a stochastic  internal- 
wave  field  on  sound  propagation.  Equations  of  the  model  are 
approximately  solved  using  a type  of  WKB(J)  expansion  in  the  non- 
dimensional  frequency  u.  Explicit  connection  is  made  between 
environmental  and  acoustic  parameters  by  assuming  that  F = Mu, 


where  M is  the  Mach  number,  is  asymptotically  order  one  with  respect 
to  u.  Other  assumptions  include  considering  only  SOFAR  rays  which 
are  regularly  perturbed  by  order  M sound-speed  and  velocity  changes. 
Our  results  modify  traditional  ray-theory  formulas  for  the  order 
one  phase,  but  not  the  amplitude,  of  a single  ray.  Both  deterministic 
and  stochastic  velocity  and  sound-speed  perturbations  may  be 
treated  by  our  model.  As  an  example  of  an  important  statistic 
of  the  received  acoustic  field,  we  present  a formula  for  the 
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expected  multipath  intensity,  assuming  stochastic  ray  phases  which 
are  both  independent  and  Gaussian. 

Interfacing  conditions  between  acoustic  and  hydrodynamic 
models  are  stated  and  motivated.  These  conditions  depend  on  the 
values  selected  for  the  variable  scales,  such  as  characteristic 
velocity  UQ  and  length  XQ.  In  this  paper,  the  choice  for  XQ  is 
ocean  depth,  but  our  techniques  may  be  applied  to  other  types  of 
problems,  for  which  different  scales  are  appropriate.  Interfacing 
trapezoids,  which  bound  permissible  regions  of  transmission  range 
Rq  and  CW  frequency  fQ,  are  constructed  and  interpreted.  As  one 
example,  the  trapezoid  in  Fig.  2 for  UQ  = 3 cm  sec-1  indicates 
that  for  Rq  = 200  km,  the  permissible  frequency  band  is  about  30 
to  560  Hz.  To  our  knowledge,  the  role  of  interfacing  in  providing 
such  constraints  has  been  overlooked  previously.  A consistent  and 
general  relation  between  the  Brunt-Vaisala  frequency  and  static 
sound-speed  distribution  is  also  developed.  Significant  simplification 
of  the  fluctuation-dependent  ray  phase  is  obtained  through  our 
consistent  approach. 

We  apply  our  model  to  the  derivation  of  formulas  for  ray  phase 
variance  in  the  presence  of  a stochastic  internal-wave  field.  The  phase 
variance  is  a fundamental  statistic  for  both  total-field  intensity  calcu- 
lations and  reciprocal  acoustic  transmissions.  To  describe  the 
internal-wave  fluctuations,  we  use  a consistent  model  previously 
published  by  the  authors.  Advantages  of  this  model  include  its 
incorporation  of  vertical  mode  structure  and  its  independence  of 
specific  internal-wave  spectral  functions  and  static-state  distri- 
butions. The  phase  variance  is  found  to  split  into  two  constituents, 
depending  on  the  velocity  and  sound-speed  fluctuations  individually. 


Approximate  evaluation  of  our  formulas  for  the  constituents 
is  feasible  for  two  classes  of  rays,  one  including  axial 
and  "near-axial"  rays  and  another  for  which  the  rays  possess 
relatively  significant  variation  in  depth  over  transmission  range. 

We  evaluate  the  constituents  for  the  second  class  using  a stationary 
phase  approximation.  Justification  for  this  procedure  is  a direct 
consequence  of  our  consistent  model  formulation.  Moreover,  the 
consistency  and  generality  of  our  analysis  enables  us  to  obtain 
many  new  results  and  predictions  for  the  variance  constituents. 

For  example,  sound-speed  (velocity)  contributions  to  the  variance 
are  apportioned  according  to  the  amount  of  internal-wave  expected 
potential  (kinetic)  energy.  The  relative  significance  of  these 
contributions  is  shown  to  be  strongly  dependent  on  a parameter  I, 
which  characterizes  the  frequency  distribution  of  internal-wave 
energy.  Further,  the  internal-wave  mode  number  spectrum  is  shown 
to  significantly  modify  the  depth  dependence  of  the  variance 
constituents.  The  degree  of  influence  depends  on  model  parameters 
and  ray  turning  depths;  in  one  case  which  can  be  compared  with 
previously-published  results,  our  predictions  range  from  60% 
greater  to  30%  less. 

We  apply  our  phase  variance  results  to  the  evaluation  of 
expected  total-field  intensity  in  an  ocean  modeled  by  a bilinear 
static  sound-speed  distribution.  Consistency  requirements  are 
used  to  determine  an  appropriate  bi-constnat  Brunt-Vaisala  frequency 
and  to  suggest  feasible  values  for  propagation  parameters.  We 
show  that  deviations  from  both  coherent  and  incoherent  intensities 
produced  by  internal  waves  can  be  very  significant,  and  that  neither 
of  these  intensities  provides  an  accurate  picture  of  the  variations 


with  Rq  in  the  presence  of  internal  waves.  The  effects  of  variations 
of  the  internal-wave  parameters  I and  j * are  examined  and  inter- 
preted. Increased  fluctuations  of  the  expected  intensity  about 
the  incoherent  are  observed  as  either  of  these  parameters  increase. 
For  the  cases  we  examined,  the  magnitudes  of  these  fluctuations 
can  be  very  dramatic,  as  high  as  21  dB  over  the  interval  of  Rg 
values  considered. 
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